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1. INTRODUCTION 

Distributed renewable energy generation is small power-generating unit clubbed with modern 
technologies and present diverse challenges while integrating with utility grid Most of the distributed 
renewable energy sources (D-RESs) shall be located near load centre on distribution grid. D-RES introduces 
negative effects on system operation. One of the negative impacts is relay coordination [1]. The integration of 
distributed renewable energy sources (D-RES) on the electrical network contains many advantages. 
However, the result of introducing these energies in the network will change the value and direction of 
energy flow as well as short-circuit current levels in different points .Several factors may lead to the effects 
of RES energies on protection devices, including the type of technology, capacitance and the position in the 
network [2]. 

The authors in [3], proposed a dual simplex method to find the setting of TDS while keeping 
constant values of Ip in modified IEEE 33-bus with installing the D-RES. In the present article, a Particle 
swarm optimization is proposed to find out the new setting for both TDS and Ip of relays with respect to 
various constraints in the presence of two D-RES. 
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2. DIRECTIONAL OVER-CURRENT RELAYS 

For optimal coordination of DOC Relays, it is necessary to determine the optimal values of TDS and 
IP of each relay in the system, this optimization of TDS and IP for minimize the operating time of the 
primary relays. The tripping time equation of DOC Relay is a non-linear equation, Depending on the fault 
current and both TDS and IP settings. Thus, the relay tripping time equation for DOC Relay is given by (1) 
[4], [5]. 
T =TDS X ae (1) 
(cerry) = 
Where, KCT is the current transformer ratio. 
The pickup value is the minimum value of the current for which the relay operates. The operating time (Tot) 
of the relay is defined by the TDS for each current value. The limits on the TDS can be presented as (2) 


TDS" < TDS; < TDs 2) 
The TDS are assumed to vary between 0.05 and 1.2. 

I” < lpi S Ip (3) 
The pickup current setting Ip are assumed to vary between 0.25 and 2.5 with step 0.25 

When a fault occurs, there must be primary protection for isolat as soon as possible to reduce the 
isolated area as much as possible, and in the case of primary relay, other relay should be protected as backup 


protection to assure the requirement of selectivity. In order to satisfy selective requirements, you must add 
the following limitation [6], [7]. 


Tact g To imäry < CTI (4) 
T — TE < CTI (5) 


Where 

Tj! is the tripping time of jth backup relay; 

Ti! is the tripping times of ith primary relay. 

The CTI is the minimum time gap in operation between the primary and its backup relay, is usually selected 
between 0.2 s and 0.5 s [4], [7], [8]. The CTI is set to a desired value of 0.3 sec. 


3. IMPLEMENTATION OF PROPOSED PSO ALGORITHM ON THE COORDINATION 
PROBLEM 
The Particle Swarm Optimization (PSO), which has gained rapid popularity as an efficient 
optimization technique, is relatively a recent heuristic introduced by Eberhart and Kennedy [9]. The various 
steps involved in the implementation of PSO to the DOCRs problem with 64 decisions variables are: 
Step1. Input the data of 33 IEEE distribution test system for calculate the power flow and current 
Step2. Define population size (50), no of iteration (=500), assume suitable values of PSO parameters 
Step3. Initial searching points and velocities are randomly generated within their limits. Pbest is set to each 
initial searching point. The best-evaluated values among individual best (P®*') are set to global best 
(gresty, 
Step4. New velocities are calculated using the equation (9). 
Step5. Evaluate the fitness values for new searching point. If evaluated values of each agent is better than 
previous Pbest then set to Pbest. If the best Pbest is better than best gbest then set to gbest. 
Step6. If the maximum iteration is reached stop the process otherwise go to step3. 
The summation of tripping times of all primary relays is known as objective function (OF). 
Conventionally, the optimal coordination problem determines two parameters, that is, the pickup current 
setting Ip and the time dial setting TDS. The OF is given by: 


OF =miny, ti (6) 


Where: OF is the objective function in sec, 





Application of PSO for optimal coordination of directional over-current Relays..., (Lazhar Bougouffa) 


190 o ISSN: 2252-8814 


ti is the operating time of the ith relay, 

N is the total DOCRs in the system. 

The two variables in the objective function are TDS and Ip. By minimizing both the TDS and Ip we can able 
to minimize the tripping time of the primary relays and satisfying the CTI. 


3.1. Goal of Study 

This work, presented the PSO algorithm to find the optimal CTI of DOC Relays. The proposed 
algorithm was evaluated using IEEE 33-bus system in presence of two D-RESs. The main contributions of 
this paper are using the optimization technique to find optimal relay settings, solve the non-linear 
coordination problem of DOC Relays with the system changed topology. 


4. CASE STUDIES 

In this paper IEEE 33-bus radial distribution feeder is considered to optimal pair of DOC Relays in 
presence of two emplacements of D-RES. The operating voltage is 12.66 kV. DOC Relays are used in the 
system to protect the feeder in case of three phase faults. The IEEE 33-bus system is presented in 
Figure 1 [10], [11]. The model is simulated in MATLAB software without and with Renewable Energy 
Sources to calculate the optimal setting of DOCRs using PSO method under system changes. The program is 
implemented for two emplacements of Renewable Energy Sources. Three phase Faults are generated on each 
bus. 





Figure 1. Power distribution system of the IEEE 33-bus system under study 


4.1. Impact of the D-RES on the short circuit level 

From a power perspective, the integration of a new D-RES into the distribution networks will 
increase the power of the network. Thus the new source will be modeled in such a way as to take only its 
current contribution in the event of a short circuit. The electrical network, of which we will determine the 
different short circuit currents without (and with) D-RES, has a radial architecture Figure 1, and we take a 
part of this network for the calculations of short circuit currents, as show in Figure 2. 
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Figure. 2. Simple distribution network with RES 


The current in the fault point will be the sum of the short circuit current without D-RES and the 
short circuit current injected into the network only by the D-RES, is given by (7) 


lec = Tora e Pe (7) 


The direct impedance seen between the bar (11) and the earth is (8) 


D-RES_ — ,D-RES ;yD—-RES 
Z a4 d-cc + jXec os Za-Line2 (8) 


Implies that (9), (10), (11). 


y? 
[D-RES = d (9) 
a = -IRES ZD RES Z 
zP RES G_ oct iXFe RES 4 74_Line2 

D-RES 2 Va 

In = a.-r eo (10) 

p ZD-RES 4 _cc+jX; RES +7 4-Line2 

D-RES Va 

In = A. ZDRS _ (11) 

e ZD-RES 4_cctjX8: RES 474 _Line2 
We find the fault current is (12), (13), and (14). 

P p 
V, Vv 
la = 4 + SIRES F D RES (12) 
Zd-cctZd-Line1tZd-Linez Z d—cct IXcc +Zd-Line2 
P p 
V, V 
L =a at pra lens | (13) 
d : = 77 D=RES 
Za-cc+Zd-Line1+Zd-Linez 2? ~FFSq_cet iXec +Zq-Line2 
P p 
V, vV 
Leä —— E + wee (14) 
, : = 7 D=RES 
Zd-cc+Zd-Line1+Zd-Line2 2? ~FFSq_cetiXec +Zq-Line2 


5. RESULTS AND DISCUSSION 

Two locations are chosen for installing the D-RES in the distribution system; to illustrate the effects 
of D-RES insertions on the setting of Relays. A minimum relay tripping time is the main goal of setting of 
relays. Also the CTI is satisfied (i.e., CTI > 0.3). Table 1 shows the optimal results of TDSs in presence of D- 
RES installation in the 33-bus distribution system. Table 2 shows the optimal results of Ips in presence of D- 




















RES. 
Table 1. Setting values of relais TDS Table 2. Setting values of relais Ip 
Setting of TDS Setting of Ip 
Relays Without D-RES With D-RES Relays Without D-RES With D-RES 
TDS, 0.5000 0.5000 Ip: 0.7300 0.7300 
TDS» 0.3500 0.3500 Ip2 0.7500 0.7500 
TDS; 0.3400 0.3400 Ip3 1.0400 1.0400 
TDS, 0.3300 0.3300 Ip4 0.8100 0.8100 
TDS; 0.3400 0.3000 Ips 0.3500 0.5000 
TDS 6 0.3500 0.3500 Ipe 0.3500 0.3500 
TDS; 0.2700 0.2700 Ip7 1.0700 1.0700 
TDS s 0.2400 0.2000 Ips 1.0100 1.2000 
TDS o 0.3100 0.3800 Ipo 0.3600 0.3600 
TDS 10 0.4600 0.5700 Ip 10 0.5200 0.5200 
TDS 11 0.4200 0.5200 Ipu 0.5600 0.5300 
TDS 2 0.3000 0.4000 Tp 12 0.5400 0.5100 
TDS 3 0.3400 0.4500 Ipis 0.5600 0.5500 
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Table 1. Setting values of relais TDS (Continued) Table 2. Setting values of relais Ip (Continued) 
Setting of TDS Setting of Ip 
Relays Without D-RES With D-RES Relays Without D-RES With D-RES 
TDS 14 0.3700 0.4100 Ip 14 0.5200 1.0500 
TDS 15 0.3600 0.4500 Ip is 0.5300 0.5000 
TDS 16 0.3000 0.3900 Ip 16 0.5000 0.5300 
TDS 17 0.3200 0.4500 Ipı7 0.8300 0.5300 
TDS ig 0.4000 0.0700 Ip is 1.2200 1.2300 
TDS 19 0.2600 0.1700 Ip io 1.0300 1.0000 
TDS 2 0.3500 0.3700 Ip 20 0.5000 0.5200 
TDS 2; 0.3200 0.3300 Ipai 0.5500 0.5100 
TDS 2 0.3300 0.0900 Ip» 0.9800 1.0300 
TDS 23 0.2900 0.0500 Ip 3 1.0200 1.0200 
TDS 24 0.2800 0.3200 Ip 24 1.0100 0.9800 
TDS 95 0.4300 0.4000 Ip 2s 0.5200 0.7500 
TDS 26 0.4400 0.3800 Ip% 0.5500 0.5200 
TDS 27 0.3200 0.2600 Ip 27 0.5700 0.5100 
TDS 2 0.3300 0.2700 Tp 28 0.5200 0.5600 
TDS 29 0.3400 0.3700 Ip 29 0.7200 0.7400 
TDS 30 0.2600 0.2800 Ip 30 1.5600 1.6000 
TDS 3; 0.3300 0.3500 Ips 1.0900 1.1300 
TDS 32 0.4200 0.4400 Ip 32 0.5100 0.5200 








The results obtained after the study from the proposed PSO Algorithm of the impact of D-RES on 
the DOC-Relays settings are displayed in Tables 1 and 2. By examining the results obtained, it can be said 
that D-RES increased the values of short circuit currents at the fault point. This increase will be reflected in 
the reduced operating time for protection before damage to the equipment that makes up the electrical 
network. 


6. CONCLUSION 

In this paper the authors presented the Particle Swarm Optimization method to explain and visualize 
the influence of D-RES on the protection settings of DOC Relays in the test system. The result of using PSO 
algorithm on an IEEE 33-bus is compared in two cases with and without D-RES. From analysis result, the 
impacts of D-RES location in the distribution system had been observed on the optimal relays settings. The 
results show the optimal settings of DOC-relays for eliminate the three fault currents. 
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